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representation of the multivariate Tutte polynomial together with the Penrose 
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1 Introduction 



A decade ago, Sokal [27] proved that if G = (V, E) is a loopless graph-*^ of maximum 
degree A, then all the roots (real or complex) of the chromatic polynomial Pcio) he 
in the disc \q\ < C'(A), where C(A) are semi-explicit constants (given by a variational 
formula) satisfying C(A) < 7.963907A.^ More generally, Sokal proved a bound on the 
zeros of the multivariate Tutte polynomial [29] (also known in statistical mechanics 
as the Potts-model partition function [25, 33, 34]) 

ZG{q,w) = Y.^'^^^U'^e (1.1) 

ACE eeA 

[here k{A) denotes the number of connected components in the subgraph (V", A)] when 
the edge weights w = {we} lie in the "complex antiferromagnetic regime" |l+t/7e| < 1: 

Theorem 1.1 [27, Corollary 5.5] Let G = {V, E) be a loopless graph equipped 
with complex edge weights w = {we}eeE satisfying |1 +We\ < 1 for all e. Then all the 
zeros of ZG{q,w) lie in the disc \q\ < KA{G,w), where 

A{G,w) = max^ltu^l (1.2) 



e3x 



and 



K = min^L: inf o;-^ V e"" L-^""^) ^ < 1 S> (1.3a) 

— mm - — J- ^ (1.3Dj 

a>o log(l-Fae-«) ^ ^ 

^ 7.963 906 075 890 002 502 ... . (1.3c) 

Moreover, we rigorously have K < 7.963907. 

Here the simpler formula (1.3b) for the constant K is due to Borgs [9, Theorem 2.1]. 

The purpose of this paper is to extend Sokal's bound by removing the condition 
that |1 -|- We| < 1 for all e. More precisely, we shall prove: 



^AU graphs in this paper arc finite and undirected; furthermore, they are allowed to contain loops 
and multiple edges unless we explicitly state otherwise. 

^More recently, Borgs [9] has provided a simpler variational characterization of the constant 
K = limA^oo C'(^)/^ ~ 7.963906 than the one given by Sokal [27, Proposition 5.4] — compare 
eqs. (1.3a) and (1.3b) below — and Fernandez and Procacci [14] have provided, in an analogous 
way, a simpler variational characterization of the constants C(A). Furthermore, Fernandez and 
Procacci [14] have improved the constants C(A) to smaller constants C*(A), for which K* = 
limA^oo C*(A)/A w 6.907652. 
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Theorem 1.2 Let G = {V, E) he a loopless graph equipped with complex edge weights 
w = {we}eeE- Then all the zeros of Zc^q^w) lie in the disc 



\q\ < /C*(*(G,to))A'(G,w;), 



(1.4) 



where 



/S.'{G,w) — max min <f Iwel, 1 — ^ — - — -\ 
^(G, 11?) = max IT max{l, |1 + Well 



e3x 



and 



ICii)) = min<jL: inf (e" - 1)"^ J] e"" V"^^ L-^""^) 



n 



n-l 



n=2 



< 1 



mm 



1<J/<1+V-V2 (1 + -0-1/2 _y) logy 



= w 



1 - W 



1 + jjj-y^ 



(1.5) 
(1.6) 



where W is the Lambert W function [11], i.e. the inverse function to x ^ xe^ 



(1.7a) 
(1.7b) 

(1.7c) 
(1.7d) 



When |1 + We| < 1 for aU e, we have A'{G, w) = A(G, w) and ^(G, w) = 1, so 

that Theorem 1.2 reduces in this case to Theorem 1.1 with an improved constant [14] 
K* ^IC*{1) = W{e/2)/[l - W{e/2)f ^ 6.907 651 697 774 449 218 .. . . This explicit 
formula for the Fernandez-Procacci [14] constant K* appears to be new. 

Let us also remark that the upper bound (1.7d) gives precisely the first two terms 
of the large--^ asymptotics of K,*{'4>): see (6.36). 

Please note that ^{G, w) involves a product over edges e3 x rather than a sum, 
and hence grows exponentially (rather than linearly) with the vertex degree whenever 
|1 + We\ > 1. The resulting exponential dependence of the bound on \q\ given in 
Theorem 1.2 is not merely an artifact of our proof, but is a genuine feature of the 
regime |1 + We| > 1.^ To see this, it suffices to note that whenever one replaces 
an edge e by k edges in parallel, the effective couplings We,eff = (1 + We)'' — 1 grow 
exponentially in k when 1 1 + j > 1 but only linearly when 1 1 + We | < 1 • For 
instance, the graph G = k':^'' (a pair of vertices connected by k parallel edges) with 
all edge weights equal has ZG{q,w) = q[q + (1 + w)^ ~ so that we must take 
|g| > 1(1 + w)'^ — 1| to avoid a root. This has roughly (but not exactly) the same 



^See also [27, Remark 2 after Corollary 5.5]. 
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dependence in w and k as the bound of Theorem 1.2. See Example 7.3 below for 
details. 

The linear growth of the bound (1.4)/(1.7) as ^{G,w) — > oo is, however, an 
artifact of allowing multiple edges. If we restrict attention to simple graphs, then 
with a little more combinatorial work we can obtain a bound that grows only like 
*(G',-i«)V2: 

Theorem 1.3 Let G — iY-iE) he a simple graph (i.e. no loops or multiple edges) 
equipped with complex edge weights w — {we}e^E- Then all the zeros of Zaiq^w) lie 
in the disc 

\q\ < Kl^i'iCwy/'AiCw), (1.8) 

where 

A(G, ly) = max min I Iwel, 1 ^ — -^ttttI (1-9) 

e3x ^ ' ei ^ 

and A = A'(G', w) /A{G, w) and 

Kl = min |l: inf (e" - 1)'^ ^ e^L-^-^) " : j^^^ < l[ (1.10a) 



oo 



.(„_,) [l + (n-l)Ar ^ 



a>Ci ^ (n — 1)\ 



n=2 



= min ^- (1.10b) 

i<y<2 {2 -y) logy 

< 5 + 2A. (1.10c) 

The constant is an increasing function of A = A'(G, w)/ A{G, w) G (0, 1], but 
the variation is fairly weak. In the complex antiferromagnetic regime |1 + We| < 1, 
where A = 1, we have — K* 6.907652, i.e. the same constant as in Theorem 1.2 
(indeed, the two theorems give the same bound in this case); while in the complex 
ferromagnetic regime |1 + We| > 1, where < A < 1, decreases as \we\ ^ oo down 
to the value = W{2e)/[2 [W{2e) - 1]^] ^ 4.892888. It is worth noting that the 
weaker version of Theorem 1.3 in which i^T^ and A{G,w) are replaced by and 
A(G, w), respectively, can be proven by a shghtly simpler argument (see the Remark 
in Section 6 after the proof of Theorem 1.3).^ 

Please observe that in the bound (1.8) we pay a price, compared to (1.4), by having 
A{G,w) in place of A'{G,w). In fact, the simple example G = K2 shows that the 
bound of Theorem 1.3 can in some cases be inferior to that of Theorem 1.2, by a factor 
of up to Kq/A 1.223222 (see Examples 7.1 and 7.2 below). But this seems to be 
the largest possible ratio of the two bounds; and in any event the ratio cannot exceed 
K*/4: 1.726913 (see the Discussion in Section 6 after the proof of Theorem 1.3). 
In most cases, however. Theorem 1.3 is a big improvement over Theorem 1.2, due 



*The first preprint version of this paper (http://arxiv.org/abs/0810.4703vl) contained only 
this weakened version of Theorem 1.3. 
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to its replacement of "^(G,w) by ^(G, it>)^/^: see, for instance, Examples 7.5 and 
7.6. Furthermore, the growth as ^'(G, lu)^/^ is essentially best possible: see again 
Examples 7.5 and 7.6. 

It is curious that the bound of Theorem 1.3 is not always better than that of 
Theorem 1.2, despite using better "ingredients" in its proof; the reasons for this will 
be discussed near the end of Section 6. It would be interesting to try to find a single 
natural bound that simultaneously improves Theorems 1.2 and 1.3. 

Please note also (see e.g. [29]) that if G is a loopless graph with multiple edges, 
then its multivariate Tutte polynomial is identical to that of the underlying simple 
graph G in which each set of parallel edges ei, . . . , in G is replaced by a single 
edge e in G with weight We = Y[i=ii^ + "^ej ^1- So one is always free to apply 
Theorems 1.2 and 1.3 to {G,w) instead of applying Theorem 1.2 to {G,w). The 
following lemma concerning the behavior of \I'(G,tu) and A\G,w) under parallel 
reduction — which will be proven at the end of Section 6 — implies that the bound 
we get by applying Theorem 1.2 to (G, w) will never be worse than the bound we 
get by applying Theorem 1.2 to (G, w). So we can find our best bound for any given 
(multi)graph G by constructing (G, w) and then taking the minimum of the bounds 
we obtain by applying (1.4) and (1.8) to {G,w). 

Lemma 1.4 Let wi, W2 & C and put ^3 = (1 + wi){l + W2) — 1- Then 

max{l, |1 + wsl} < max{l, |1 + Will max{l, |1 + W2I} (l-H) 

and 

fi h-'-'sl 1 fi I b-'-'il 1 fi I I'"''2| 1 

mm <^ \w-i\, —— — - \ < mm <^ Wi , — }■ + mm <^ W2 , -r— r > ■ (1.12) 

[ II + W3IJ [ |l + Wi|J [ II + W2IJ 

Sokal's proof of Theorem 1.1 involved the following steps: 

1. Write the multivariate Tutte polynomial ZG{q,w) as the partition function of 
a polymer gas with weights depending on q and w (this is easy: see Section 2 
below). 

2. Invoke the Kotecky-Preiss [20] condition for the nonvanishing of the partition 
function of a polymer gas. 

3. Control the polymer weights by bounding sums over connected subgraphs by 
sums over trees, using the Penrose inequality [24]. This step required |l + We| < 
1. 

4. Bound the total weight of n-vertex trees (or more generally, of connected sub- 
graphs with m edges) in G that contain a specified vertex x & V. 

5. Put everything together to prove that Zoiq^w) ^ whenever q lies outside a 
specified disc. 
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Here we follow the same outline, but modify step 3 so as to allow arbitrary com- 
plex weights Wg. In addition, in step 2 we replace the Kotecky-Preiss condition 
by the more powerful Gruber-Kunz-Fernandez-Procacci [16, 13] condition, thereby 
slightly improving the numerical constant along the lines of the work of Fernandez 
and Procacci [14] for chromatic polynomials. Finally, Theorem 1.3 needs a slightly 
strengthened version of the bound in step 4. 

The plan of this paper is to treat each of these five steps in successive sections. 
Thus, in Section 2 we recall how the multivariate Tutte polynomial Zciq, w) can be 
written as the partition function of a polymer gas. In Section 3 we recall the Kotecky- 
Prciss and GrTibcr-Kunz-Fernandez-Procacci conditions for the nonvanishing of the 
partition function of a polymer gas. In Section 4 we recall the Penrose identity [24] and 
show how to use it to bound the polymer weights without assuming that |1 + ti;e| < 1; 
this is our main new contribution. In Section 5 we recall the sharp bound [27, 17] in 
terms of maximum weighted degree A(G, w) on the total weight of connected m-edge 
subgraphs in G that contain a specified vertex x; we also strengthen it slightly by 
taking specific account of the edges incident on x. In Section 6 we put everything 
together to prove Theorems 1.2 and 1.3; we also prove Lemma 1.4. Finally, in Section 7 
we examine some examples that shed light on the extent to which Theorems 1.2 and 
1.3 are sharp or non-sharp. 

2 Polymer-gas representation of Zciq^w) 

In this section we recall how to rewrite the multivariate Tutte polynomial Zciq, w) 
as the partition function of a polymer gas living on the vertex set of G. This easy 
result is due to Sokal and Kupiainen [27, Proposition 2.1]. 

First, some notation: If if = (V, E) is a graph equipped with edge weights 
w — {we}eeEj we denote by Ch{w) the generating polynomial of connected span- 
ning subgraphs of H, i.e. 



Note that Ch{w) = if if is disconnected. 

If G = {V, E) is a graph and S C V, wc denote by G[S] the induced subgraph of 
G on S, i.e. G[S] is the graph whose vertex set is S and whose edges consist of all 
the edges of G both of whose endpoints lie in S. 

Proposition 2.1 (polymer representation of the multivariate Tutte polynomial) 

Let G — {V, E) he a loopless graph equipped with edge weights w = {we}eeE- Then 






(2.1) 



ACE 
(V,j4) connected 



oo 



n 





(2.2) 



N=0 {Si,...,5jv} i=l 
disjoint 
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where the sum runs over unordered collections l^*!, . . . , ^'jv} of disjoint nonempty 
subsets ofV, and the weights ^{S) are given by 

(q-(\sh^)CG\s](w) if\S\>2 
^(S) = \ ''^^^^ ' ' - (2.3) 

\0 if\S\^l 

[The N = term in the sum (2.2) is understood to contribute 1.] 

The identity (2.2) thus represents q~\^^ZQ{q, w) as the partition function of a gas 
of nonoverlapping "polymers" Uving on V , with weights (2.3). Here a "polymer" is, 
in principle, any nonempty subset 5 C \/; but since the weight S,{S) vanishes for sets 
of cardinality 1, we can equivalently restrict our polymers to be subsets of cardinality 
at least 2. Likewise, the weight ^{S) vanishes whenever the induced subgraph G[S] 
is disconnected; so we can, if we wish, restrict our polymers to be sets S for which 
G[S\ is connected. 

Proof of Proposition 2.1. Starting from the definition (1.1) of Zg(5, i^), let us 
separate the terms in the sum according to the number k of connected components 
[i.e. k{A) = k] and according to the partition {Si, . . . , Sk} of V that is induced by the 
vertex sets of those connected components; we will then sum over all ways of choosing 
edges within those vertex sets so as to connect those vertices. We thus have 

k 

ZG{q,w) = q^^'^Y. E n^"^'^''"'^^G[5d(«^), (2.4) 

k>l {Si,...,Sk} i=l 
V=\i)Si 

where the sum runs over all unordered partitions {^i, . . . , Sfc} oi V into nonempty 
subsets, and we have used \V\ — Yli=i Note now that any set Si of cardinality 1 
gets weight q^^^^^^^^^CG[s^]{^) = 1 (here we have used the fact that G is loopless). 
So let us define {SJ, . . . , S^} to be the subcollection of {S*!, . . . , iS^.} consisting of 
the sets of cardinality > 2; and let us note that there is a one-to-one correspondence 
between unordered partitions {Si, . . . , Sk} oiV into nonempty subsets and unordered 
collections {S'l, . . . , S'j^} of disjoint subsets of V of cardinahty at least 2 (which need 
not cover all of V: indeed, the points not covered correspond to the singleton sets Si 
in the original partition). Passing to {S[, . . . , S'j^f} and dropping the primes, we have 
(2.2)7(2.3). □ 



3 Sufficient condition for the nonvanishing of a 
polymer-gas partition function 

Let y be a finite set, and let {p{S)}0^scv be a collection of complex weights 
associated to the nonempty subsets of V. Consider now a gas of nonoverlapping 
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"polymers" living on V, with weights p{S): the partition function of such a polymer 
gas is, by definition, 

oo n 

^ = E E Upi^^h (3.1) 

N=0 {Su-,Sn} i=l 
disjoint 

where the sum runs over unordered collections {Si, . . . ,Sm} of disjoint nonempty 
subsets of V, and the N = term in (3.1) is understood to contribute 1. The following 
proposition — essentially proven almost four decades ago by Gruber and Kunz [16, 
Section 4, cf. eq. (33)] but largely forgotten, and then rediscovered very recently by 
Fernandez and Procacci [13, eq. (3.17)] with a new proof — gives a sufficient condition 
for the nonvanishing of a polymer-gas partition function: 

Proposition 3.1 (Gruber Kunz— Fernandez— Procacci condition) Let V be a 

finite set, and let {p{S)} 0^scv be complex weights associated to the nonempty subsets 
ofV. Suppose that there exists a number a > such that 

sup Ve"l^' |p(5)| < e"-l. (3.2) 

Then 

oo n 

^ ^ E E Up(^^) + 0- (3-3) 

iV=0 {5i,...,5jv} i=l 
disjoint 

See also [6] for an extremely simple proof of Proposition 3.1 by induction on V . 

In the slightly less powerful Kotecky-Preiss [20] condition, the term — 1 on the 
right-hand side of (3.2) is replaced by a. 

Remark. Suppose that (as happens in all nontrivial cases) there exists a set S 
with IS"! > 2 and p{S) ^ 0. Then the hypothesis that there exists o; > such that 
(3.2) holds can be rewritten as 




since in this case the infimum on the left-hand side of (3.4) will always be attained 
at some a > 0.^ We will use the Gruber-Kunz-Fernandez-Procacci condition in the 
form (3.4). 

^If there exists a set S with |5| > 2 and p{S) ^ 0, then the function f(a) being minimized on 
the left-hand side of (3.4) is a continuous function that tends to +oo as a J, and as a t hence 
its minimum is attained. 

There is one exceptional case in which (3.4) holds but there does not exist a > such that (3.2) 
holds: namely, if p{S) = whenever \S\ > 2 and in addition we have max|p({a;})| = 1. Indeed, if 

x&V 

p(S') = for \S\ > 2, we have S = H [^ + P({^})]' which vanishes when at least one p{{x}) equals 
— 1; so (3.4) fails (barely) to imply S ^ in this case. 
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4 A bound on Cij(iy) via the Penrose identity 



In this section we recall the Penrose identity [24] and show how it can be used to 
bound a sum over connected subgraphs by a sum over trees even in the absence of 

the hypothesis \1 + Wf.\ < 1. 

Let H = (V, E) be a graph. Recall that Ch{w) denotes the generating polynomial 
of connected spanning subgraphs of H: 

ch{w) ^ j2 n^- (4.1) 

ACE eeA 
(V,A) connected 

We denote by Th{w) the generating polynomial of spanning trees in H: 

Th{w) = 5^ (4.2) 

ACE eeA 
(V,A) tree 

Let C (resp. T) be the set of subsets ACE such that (V, A) is connected (resp. 
is a tree). Clearly C is an increasing family of subsets of E with respect to set- 
theoretic inclusion, and the minimal elements of C are precisely those of T (i.e. the 
spanning trees). It is a nontrivial combinatorial fact — apparently first discovered 
by Penrose [24] — that the (anti-)complex C is partitionable: that is, there exists a 
map R: r ^ C such that R(r) D T for all T e T and C = ISreriT^ ^(^)] (disjoint 
union), where [Ei,E2] denotes the Boolean interval {A: Ei C. A C. E2}. In fact, 
many alternative choices of R are available®, and most of our arguments will not 
depend on any specific choice of R. An immediate consequence of the existence of R 
is the following simple but fundamental identity: 

Proposition 4.1 (Penrose identity [24]) Let R: T ^ C he any map such that 
R(T) D T for allT E T andC is the disjoint union of the Boolean intervals [T, R(r)] 
over T E T. Then 

ch{w) = e n^^ e n (4-3a) 

TCE eeT TCACR(T) eeA\r 
(V,T) tree 

= E n^^ n (i+^e). (4.3b) 

TCE eeT eeR(T)\T 
(V,T) tree 

If |1 + We| < 1 for all e, then it is obvious that we can take absolute values 
everywhere in (4.3b) and drop the factors |1 + We|, yielding: 



^See e.g. [24], [7, Sections 7.2 and 7.3], [36, Section 8.3], [15, Sections 2 and 6], [5, Proposition 
13.7 et seq.], [27, Proposition 4.1] and [26, Lemma 2.2]. 
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Proposition 4.2 (Penrose inequality [24]) Let H = (V, E) be a graph equipped 
with complex edge weights w = {we}eeE satisfying |1 + We\ < 1 for all e. Then 

\Ch{w)\ < Th{\w\) . (4.4) 

Remark. By using a specific choice of tlie map R (namely, that of Penrose [24]), 
Fernandez and Procacci [13] have recently shown how to improve Proposition 4.2 when 
We e {—1,0} for all e; and this improvement plays a key role in their proof of the 
Gruber-Kunz-Fernandcz-Procacci condition (Proposition 3.1) for polymer gases with 
hard-core repulsive interactions. See also Fernandez et al. [12] for a generalization 
to —1 < We < 0, which leads to an improved convergence criterion for the Mayer 
expansion in lattice gases with soft repulsive interactions. □ 

Let us now show what can be done without the hypothesis |1 + We| < 1. 

Proposition 4.3 (extended Penrose inequality) Let H — (V, E) be a graph equipped 

with complex edge weights w — {we\eeB- Then 



\Ch{w)\ < TH{w')Y[ui^.^{l,\l + We\] (4.5a) 

eGE 

< Th{w')^>{H,w)\''\I\ (4.5b) 

where 

w'e = min||we|, iJ^'^ 1 1 (4.6) 

and 

'^{H,w) = maxTTmaxjl, |1 + We|} . (4.7) 



Proof. Let us write E_,_ = {e G E: |1 + We| > 1}- Taking absolute values in the 
Penrose identity (4.3b), we obviously have 

\ch{w)\ < Yi ni^^i n ii+^ei- (4.8) 

TCE eer eeR(r)\T 
(V,r) tree 

Now, in each summand on the right-hand side of (4.8), multiply and divide by jH-Wel 
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for each e e T n E+; we have 

\ch{w)\ < y: iih n jrhri n ii+-ei n \^+ 

TCE eeT eernE+ I ^' eeTnE+ eeR(r)\r 

(V,T)tree (4.9a) 

- E n--{i-ei,^^| n (4.9b) 

TCE eeT I ^ ee(TnE+)U(R(T)\T) 

(V,T) tree 

TCE eeT I ^ eeE 

(V,T) tree 

This proves (4.5a). Then (4.5b) is a trivial consequence. □ 

If we assume that the graph H is simple (i.e. has no loops or multiple edges), then 
we can get a shghtly better bound: 

Proposition 4.4 (extended Penrose inequality for simple graphs) Let H = 

(V, E) be a simple graph (i.e. no loops or multiple edges) equipped with complex edge 
weights w — {wgleeE- Then, for any vertex x &\/, we have 

\Ch{w)\ < Th{w") W max{l,|l + we|} (4.10a) 

eeE\E(a;) 

j| maxjl, |1 + j^/'' 

eeE(a;) 

< TH{w)<^{H,wf'^\-^^l'' (4.10c) 

< r^^(|Mj|)*(ii-,'u;)(l^l-i)/2 (4_;L0d) 

where E(x) denotes the set of edges incident on x, 

I \we\ if e G Efx) 

< - { (4.11) 
[minjlwel, |we|/|l + Well otherwise 

{minjlwel, |'?i'e|/|l + 'W^el^''^} if e e E(x) 
(4.12) 
min{|we|, |we|/|l + Well otherwise 



and 



-^{H^w) = maxTTmax{l, |1 + We|} . (4.13) 

e3y 
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It is worth observing that (4.10a) is a genuine improvement of (4.5a), because 
the product fl^^^^^^ max{l, |1 + We\} more than compensates the factors w"/w'^ = 
max{l, ll+Wel} for the subset of edges in E(x) that happen to he in any given spanning 
tree of H [that is, we have Th{w") < Th{w') JJ^g^^^^ max{l, |1 + Wgl}]. Similarly, 

(4.10c) is always an improvement of (4.5b), since Th{w) < Th{w') 'if{H,wy^^. 

The main change in Proposition 4.4 compared to Proposition 4.3 is that the power 
of "${11, w) is reduced from |V|/2 to (|V| — l)/2. It is perhaps rather surprising 
that such a small modification in this intermediate bound results in the significant 
improvement obtained in going from Theorem 1.2 to Theorem 1.3 — namely, reducing 
the growth from ^(G, w) to its square root — but that is how things turn out. An 
explanation will be given in Section 6, after the proof of Theorem 1.3. 

We do pay a price for this improvement: namely, the argument of Th in (4.10c) 
is w rather than w'. This is what causes the bound of Theorem 1.3 to depend on 
A{G,w) rather than the somewhat smaller A'(G,w). We do not know whether this 
can be improved. 

The proof of Proposition 4.4 — unlike all the preceding results in this section — 
depends on a specific choice of the map R, namely the one used by Penrose in his 
original paper [24]. Let us briefly recall Penrose's construction (see [13, 12] for more 
details). We assume that H — (V, E) is a simple graph, and we choose (arbitrarily) 
an ordering of the vertex set V by numbering the vertices 1, 2, . . . , n (where n — |V|). 
We consider the vertex 1 to be the root, and denote it by r. If T C E is the edge 
set of a spanning tree in H [that is, (V, T) is a tree], then for each a; G V we denote 
by distT(a;) the graph-theoretic distance in the tree (V, T) from the root r to the 
vertex x. Given T, the vertex set V is thus partitioned into "generations" , defined as 
the sets of vertices at a given distance from the root r. 

The Penrose map R: T ^^ ^(T) is then defined as follows. For any tree T C E, 
the edge set R(r) D T is obtained from T by adjoining all edges e e E that either 

(a) connect two vertices in the same generation [i.e. at equal distance from the 
root r in the tree (V, T) — note that no such edge can belong to T] , or 

(b) connect a vertex a; to a vertex x' in the preceding generation [i.e. with dist-i (x') = 
distr(a;) — 1] that is higher-numbered than the parent of x [here the parent of 
X is the unique vertex y with distr(y) = distT(a;) — 1 such that xy & T]. 

It can be shown [24, 13, 12] that R is indeed a partitioning map in the sense that C is 
the disjoint union of Boolean intervals [T, R(T')]. Furthermore, it follows immediately 
from this construction that R(T) \ T cannot contain any edge incident on the root r; 
that is, R(T) \ T C E\ E(r).^ 

''We remark that this would no longer be the case in a generalization to the Penrose construction 

to non-simple graphs. In such a generalization, wc would also order the edges connecting each pair 
of vertices, and wc would add to the definition of R(T') a third case: 

(c) connect a vertex x to its parent y by any edge that is higher-numbered than the edge con- 
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With the fact R(r) \ T C E \ E(r) in hand, the proof of Proposition 4.4 is almost 
immediate: 

Proof of Proposition 4.4. We choose an ordering in which the vertex a; G V 
appearing in the statement of the proposition is numbered 1, i.e. r = x. Now follow 
the proof of Proposition 4.3 through (4.9b). We have (T n E+) U (R(T) \ T) C 
(E \ E(a;)) U (T fl E+ n E(a;)). The terms in E \ E(x) can be bounded above by 
neeE\E(a;) m^xfl; |l + We\} analogously to what is done in (4.9c), and the terms in 
T n E_|_ n E(x) convert the argument of Th from w'^ — mm{\we\, \we\ / \1 + We\} to w" . 
This proves (4.10a). 



Then 



n n max{l,|l + «;e|}V2 
max-^1, l + Wef = 7 — \ tttit; — (4.14a) 



eeE\E(a:) 



n max{l, |l + We|}i/2 

^(H ti,)(|v|-i)/2 
< ^[-n,w) 

- n max{l, |l + We|}V2 ^ ^ 

eeB{x) 

since the numerator of (4.14a) counts every edge in E \ E(a;) twice and every edge in 
E(a;) once; so (4.10b) follows from (4.10a). Moreover, by definitions (4.11), (4.12) and 
(4.2) we have 

< Th{w) , (4.15) 



THiw") 



n max{l, |1 + We|}]i/2 



SO (4.10c) follows from (4.10b). Finally, (4.10d) is a trivial weakening of (4.10c). □ 



5 Bounds on connected m-edge subgraphs con- 
taining a specified vertex 

Let G = {V, E) be a graph equipped with edge weights w = {we}eeE- Let us 
define the weighted sum over connected subgraphs G' — {V, E') C G that contain 



necting a; to y in T. 

We would then no longer be able to guarantee that R(T) \ T contains no edges incident on the 
root r; rather, we could assert only that R(T) \ T cannot contain any edge incident on the root r 
that is the lowest-numbered among its set of parallel edges. 
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the specified vertex x and have exactly m edges: 



G'={V',E')CG eeE' 
G' connected 
V'3x 
\E'\=m 

We will abbreviate Cm{x;G,w) to c^(x) when it is obvious which weighted graph 
{G,w) we are referring to. We then have the following sharp bound on Cm{x) in 
terms of the maximum weighted degree 

A{G,w) = max^|we|. (5.2) 

Proposition 5.1 [27, Proposition 4.5] Let G = {V,E) be a loopless graph 
equipped with edge weights w — {we}eGE- Then for any x &V, we have 

cm{x) < ^^^^^ A(G,wr (5.3a) 

m! 

< e"'A{G,wr. (5.3b) 
[The e in (5.3b) denotes, of course, the base of the natural logarithms!] 

See also [17, Section 6] for two alternate proofs of a slight generalization of Propo- 
sition 5.1, as well as examples showing that it is (in a certain sense) sharp. 

We shall see in the next section that Proposition 5.1, combined with the results of 
Sections 2, 3 and 4, is sufficient to prove Theorem 1.2 (as well as a shghtly weakened 
version of Theorem 1.3). In order to prove the full Theorem 1.3, we shall need a slight 
improvement of Proposition 5.1 in which the edges incident on the "root" vertex x 
are treated differently from the remaining edges. Let us write E{x] G) for the set of 
all edges in G that are incident on x, and define the weighted degree at x by 

d{x;G,w) = ^ \we\ . (5.4) 

eeE{x;G) 

We will abbreviate E{x; G) to E{x) when it is obvious which graph G we are referring 
to. We then have: 

Proposition 5.2 Let G — (V^, E) be a graph equipped with edge weights w — {wejeeE 
and let x eV. Put d — d{x; G, w) and D — A(G — x, w\ g-x)- Then 

d{d + mDr-' 

Cm{x) < ■ . (5.5) 

m! 



14 



Since d and D are both bounded above by A(G, lo), this is a genuine improvement 
of Proposition 5.1. 

In the proof of Proposition 5.2 it will be convenient to employ the quantities 

, \ Kim -\- / m\ for m > 1 , ^, 

C[m,K,) = < (5.6) 
11 for m = 

defined for integer m > and real k. Then (5.5) can be rewritten in the form 

Cm(x) < C{m,d/D)D'^ . (5.7) 

One way of proving Proposition 5.2 is to use Lemma 5.3 below together with 
Proposition 5.1 applied to G — x. But it is instructive to give here an ah initio proof, 
by induction on m, that is similar to the first proof of Proposition 5.1 given in [17, 
Section 6]. It relies on the following properties of C(m, n): 

(a) For each integer m > 0, C(m, k) is a polynomial of degree m in k, with nonneg- 
ative coefficients. In particular, C{m,K) is an increasing function of k for real 
K > 0. 

(b) Generating function: If C{z) solves the equation 

C(^) = e"^(") , (5.8) 

then 

oo 

C{zY = Y^C{m,K)z'^ (5.9) 

m=0 

for all real /t; this follows from the Lagrange inversion formula. Moreover, the 
series (5.9) is absolutely convergent for \z\ < 1/e and satisfies C(l/e) = e. 



(c) For integer A; > 1, 



C{m,k) = \{C{mi,l). (5.10) 



mi,...,mfc>0 i=l 

miH hntfi. =m 



This is an immediate consequence of (5.9). 
(d) For all real k and z, 



m f 



C(m,«) = ^ C(m-/,«-z + /). (5.11) 
/=o ■ 



See [17, eq. (6.7)]. 



15 



For any subset F (1 E, we use the notation w{F) = JleeF^e- We will need two 
further results from [17]: 

Lemma 5.3 [17, Section 6, Facts 1 and 2] Let G — (V, E) be a graph equipped 

with nonnegative real edge weights w = {we}eeE CLnd let x eV . For each F C E{x), 
let = {x{ , , . . . , x^(^p)} be a labeling of the vertices of V — x that are incident 
with edges in F . Then, for all m > 1, 

m 

Cm{x]G,w)< ^ w{F) Y\^Cmi{xf]G - x,w\g-x) ■ (5.12) 

0^FCE(x) mi,...,mj(p)>Q i=l 
miH \-mj(^F)=m-\F\ 



Please observe that j{F) < \F\] and if the graph G is simple, then j{F) = \F\. Please 
note also that Fact 2 as stated in [17] writes Cmi{xf ;G,w) in place of Cmi{xf ;G — 
X, w\g-x) in (5.3), but the improved version stated here is an immediate consequence 
of the proof given there. 

Lemma 5.4 [17, Lemma 6.2] Let S be a set in which each element e E S is given 
a nonnegative real weight Wg. Then, for each integer f >0, we have 

En^e < ^ fe^e) . (5.13) 

FCS eeF ■' ' \eeS / 

\F\=f 



Proof of Proposition 5.2. We will prove (5.5)/(5.7) by induction on m. The 
statement holds trivially when m = 0, so let us assume that m > 1. By Lemma 5.3, 

m 

Cm{x) < ^ \w{F)\ ^ Yl^rni{Xi;G - X,w\g-x) 

0^FCE{x) mi,...,mj(p)>Q i=l 
miH \-mj(p)=m-\F\ 

0^FCE(x) mi,...,mj(F)>Q i=l 
miH \-mj(^F)=m—\F\ 

= \w{F)\C{m-\Flj{F))D^-\^\ 

0ytFCE{x) 

< Yl \w{F)\C{m-\F\,\F\) D"'-\''\ (5.14) 

0ytFCE{x) 
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where the second hne used the induction hypothesis (5.5) apphed to the graph G — x 
(note that rui < m)^ and the fact that d{v; G — x, w\g-x) < D ioi all v & V — x, the 
third hne used the identity (5.10), and the last line used j{F) < \F\ and the fact that 
C(m, k) is an increasing function of k. Using Lemma 5.4, we have 

Ux) < D-Y.^^^C{m-fJ) 
f=o ^' 

= G{m,d/D)D'^ , (5.15) 

where the second line used C(m, 0) = for m > 1, and the last line used identity 
(5.11) with K = z = d/D. This proves (5.7). □ 



6 Proof of Theorems 1.2 and 1.3 and Lemma 1.4 

We can now put together the results of the preceding sections to prove Theo- 
rems 1.2 and 1.3. At the end of this section we will also prove Lemma 1.4. 

Proof of Theorem 1.2. We want to show that Zaiq^w) ^ whenever \q\ > 
/C*(^(G, w)) A'{G, w). We wiU do this by verifying the condition (3.4) for the poly- 
mer weights (2.3), which we recall are 

^{S) = g-(l^l-^)CG[5](t«^) for 1-51 > 2 . (6.1) 

By Proposition 4.3 [bound (4.5b)] applied to H = G[S], we have 

\Cg[s]{w)\ < TG[s]{w')^{G[Siwf\/' (6.2a) 

< TGis]iw')^iG,wp/'' (6.2b) 

where w' is defined by (4.6). Then by Proposition 5.1 apphed to the weights w', and 
observing that the n- vertex trees are a subset of the connected graphs with m = n — 1 
edges, we have for any vertex x &V 

Y.'^oidW) < ^A(G,wr-' (6.3a) 

\S\=n 

„n-l 

= —^A'{G,w)''-K (6.3b) 



^Actually, we used here only the weakening (5.3a) of the inductive hypothesis (5.5), apphed to 
G-x. 
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Therefore, the condition (3.4) for the weights (2.3)/ (6.1) is verified as soon as 



inf(e°-l)-'^e""|gr("-^)— pA'(G',M;)"-^*(G',w;)"/2 < 1. (6.4) 

n=2 

If we set L = \q\/A'{G,w), this is precisely the inequahty contained in the right- 
hand side of (1.7a). So Zciq^w) ^ whenever L > )C* {'^ {G , w)) , i.e. whenever 
\q\ >lC*{<i!{G,w))A'{G,w). 

The equivalence with (1.7b,c) and the inequality (1.7d) are proven in Lemma 6.1 
and Corollary 6.2 below: see (6.13a-c). □ 



Lemma 6.1 For A > and P > 0, define the function 

FM = min|L: inf (e" - 1)-^ V e"" L'^^''^ + ~ < (3 

I a>o^ (n-1)! ~ 



(6.5) 



Then: 

(a) Fx{(3) is an increasing function of X and a decreasing function of j3. 
(h) (5F\{(5) is an increasing function of both A and (5. 

(c) Fx[ii/\)/\ is a decreasing function of both \ andii{> 0). In particular, Fx{P)/X 
is a decreasing function of both A and (5. 

(d) logFx^P) is a convex function o/log/?. 

(e) We have 

Fx((3) = min . (6.6) 

-^^^^ i<y<i+p {1 + (3 -y)logy ^ ' 

(f) For A = 0, 1 we have 

m) = ^W^((l + /3)e)/r((l + /3)e) - If (6.7) 



Fi(/3) ^ I3W 



e 



l + (3 



1 - W 



e 



(6.8) 



where W is the Lambert W function [11], i.e. the inverse function to x ^ xe^ 
(g) For Q < \ < X' we have 
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(h) For < X < 1 we have 



FxW) < 4/3-1 + (1 + 2A) 



Corollary 6.2 The function 



K:{^) = min<jL: inf (e" - 1)-^ ^ e"" V'"/^ L" 



(n-l) 



n 



n=2 



ni 



< 1 



defined for ip > satisfies 



and hence 



mm 



i<j/<i+V-V2 (1 + -0-1/^-?/) logy 



< + 3V^i/^ . 



1-W 



1 + v'-l/2 



(6.10) 

(6.11) 

(6.12) 

(6.13a) 

(6.13b) 
(6.13c) 



Only parts (e,f,h) of Lemma 6.1 will actually be used in the proofs of Theorems 1.2 
and 1.3. But we think it worthwhile to collect here some additional properties of the 
function Fx{f3): some of these will be invoked in the Discussion after the proof of 
Theorem 1.3, while others may end up playing a role in future work. 

Proof of Lemma 6.1. (a) It is immediate from the definition (6.5) that Fx{(3) is 
increasing in A and decreasing in (3. 

(b) The change of variables L' = j3L in (6.5) shows that 



(5Fx{(5) - min 



[l' : inf (e" - 1)"^ V e"" (L')-^"-^) ^ ^J" ' 



< 1 

(6.14) 



is increasing in both A and j3. 

(c) The change of variables L" = L/A in (6.5) shows that 



A 



mm 



L": inf (e"-l)-i;f^e""(L")-("-i) ' ^J^^" ' < 



n=2 



{n-l)\ 



is decreasing in both A and ^. 

(d) Suppose that we have triplets {ai,Li,Pi) satisfying 



J2 e"'" L 



a.:nr-(n-l) [1 + (n - 1)A] 



n-2 



n=2 



(n- 1)! 



< A (e"^ - 1; 



(6.15) 



(6.16) 
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for i — 1,2. Now let k G [0, 1] and define 



a — Kai + (1 — K)a2 (6.17a) 
L = L^Ll-" (6.17b) 
/9 = P^Pl-' (6.17c) 



Then Holder's inequality with p = l/n and q = 1/(1 — k) yields 
J2 e"- Z-("-i) + ["^ ^ ]\^^^~^ < P (e"i - 1)'' (e"^ - 1 



n=2 



(6.18) 



And since the function a i-^ log(e° — 1) is concave on (0, oo), we have (e"i — l)'^(e°2 — 
l)i-« < - 1. This proves (d). 

(e) The proof that (6.5) is equivalent to (6.6) will be modelled on an argument of 
Borgs [9, eq. (4.22) ff.], who proved a related result. 

Note first that c ^ ce~'^ maps the interval [0, 1] strictly monotonically onto the 
interval [0, 1/e]; and recall [31, p. 28] that its inverse map is the tree function 



^(^) = E^^"' (6-19) 



n=l 



which is convergent and monotonically increasing for < x < 1/e and satisfies 
T{ce~'^) = c for < c < 1. Moreover, it is well known (see e.g. [11, eq. (2.36)]) that 
for all real k > one has [cf. (5.9)] 

^ f^<-^+jr-\^ (6.20) 

^ m=0 

(this is an easy consequence of the Lagrange inversion formula). Writing for conve- 
nience U{z) = T{z)/z, we therefore have 

n=l ^ ' 

for all real A > 0. 

The inequality on the right-hand side of (6.5) is then equivalent to the statement 
that Xe" / L < 1/e (otherwise the sum would be divergent) and 

e°C/(Ae"/L)i/^ - e" < /3(e° - 1) . (6.22) 

Eliminating L in favor of a new variable c defined by Ae"/L = ce^'^ with < c < 1, 
and using the fact that t/(ce^^) = e^, we see that the inequality on the right-hand 
side of (6.5) is equivalent to 

c < minjl, Alog[l + /?(l-e-")]| . (6.23) 
20 



Since L = Ae"/(ce '^), and ce increases monotonically with c for < c < 1, we 
deduce that (6.23) is equivalent to 

r e"[l + /?(l-e-°)f ^f^n_,-„)<,VA_i 
L > log[l + /?(l-e-")] ^'f'^' ' ^g_24) 

[Ae"+^ if /3(1 - e-°^) > eV^ - 1 

Changing variables from Q;toy=l + /3(l — e~"), we can rewrite this as 

if 1 < 7/ < min(e^/^, 1 + /?) 



[l + p-y)logy 
A/3e 



(6.25) 



if e^/^ < y < 1 + /3 



Now we can optimize over y: the minimum will always be found in the interval 
1 < y < e^^^, so we have 

Py^ (3y^ 
FxiP) = min — = min — , (6.26) 

l<2/<mm(ei/A,i+/3) (1 + /J - log y l<2/<l+/3 {1 + P - y)\ogy 

where the final equality results from the fact that |/'^/[(l + P — y) logy] is increasing 
for e^^^ < y < 1 + p. This proves the equivalence of (6.5) with (6.6) for A > 0; and 
the case A = follows by taking limits (or by an easy direct proof). 

(f) For A = 0, simple calculus shows that the minimum in (6.6) is attained at 
y = {1 + P)/W{{1 + P)e), so that Fo(/3) is given by (6.7). Likewise, for A = 1, simple 
calculus shows that the minimum in (6.6) is attained a,t y — {1 + P) W{e/{1 + P)), so 
that Fi{P) is given by (6.8). 

(g) To prove the comparison inequality (6.9), it suffices to observe that whenever 
< A < A' and n > 2 we have 



l + (n-l)A7 - V1 + 2A' 

(just consider n — 2 and n > 3 separately). Inserting this into the definition (6.5) 
yields (6.9). 

(h) To prove the upper bound (6.10), write y — 1 + x in. (6.6) and use the inequal- 
ities 

^ < - + - (6.28) 



log(l + x) X 2 

(1 + x)^ < 1 + Ax (6.29) 
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which are vahd for all x > and < A < 1. Therefore. 



(1 + (3 -y)\ogy (3-x 

The latter function is minimized at x = (-2+a/4 + (2 + 4A)/3 + 2A/32)/[l+(2+/3)A] e 
(0, /?), with minimum value 

^ + A + I + |V(l + /3/2)(l + A/3). (6.31) 

This, in turn, is bounded above by 4,/?"^ + (1 + 2A) on the entire interval < /? < oo.^*^ 
[Alternatively, it suffices to make this proof for A = 1 and then invoke (6.9) to deduce 
the resuh for < A < 1.1 □ 



Proof of Corollary 6.2. An easy calculation proves the formula (6.12) relating 
lC*{ij) to Fi(/3), so that (6.13a-c) follow from (6.6), (6.8) and (6.10). □ 

Remarks. 1. The proof of Lemma 6.1(e) becomes a bit simpler for (3 < e^/^ — 1, 
since we then always have — e~") < e^^'^ — 1 and hence we need not worry about 
the second case in (6.24) and (6.25). This simplification applies in particular when 
A < 1 and /5 < 1, which covers what is needed in the proofs of both Theorem 1.2 
(A = 1, /3 = tj;-^/^ < 1) and Theorem 1.3 (0 < A < 1, /3 = 1). 

2. We can compute the small-/9 asymptotics of Fx{P) by expanding (6.6) in powers 
oi y — 1: the minimum is located at 

^ 2^ 16 ^ 192 ^ 3072 ^ ^ ^ 

and we have 

Fm - 4r' + (1 + 2A) - '-^^^^P + ll+26A-m^-8A3 ^. ^ 



For A = 0, 1 an alternate method is to expand (6.7)/(6.8): we obtain 



(6.33) 



Fo(/3) = 4(3-' + 1-1(3+ —p^ - -^(3^ + -^(3^ - ... (6.34) 

"^^^ ^ 48^ 192^ 15360^ 36864^ ^ ^ 

FA(3) = 4/3-1 + 3 - — /5 + —13^ - -^p^ + _ _ _ _ (g_35N 

^ 48^ 192^ 15360^ 184320^ ^ ^ 



^Proof of (6.28): Write t = log(l + x) > 0; then (6.28) states that 1/t < l/(e* - 1) + 1/2. 
This is trivially true for t > 2; and for < f < 2 it is equivalent to e* — 1 < t/{l — t/2), which is 
obvious from the Taylor series. 

^°Proof: We have 

v/(l + ci/3)(l + C2/3) < 1 + 

for all ci, C2, /3 > 0, as is easily seen by squaring both sides and using the arithmetic-geometric- mean 
inequality y^cici < (ci + C2)/2. 
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Therefore, the large- -0 asymptotics of /C*('0) = ^/^) is 

= 4^ + 3^1/2 _ 1 + _ + _____ (6.36) 

^ ^ 48 192^ 15360 184320^ ^ ^ 

3. We conjecture (based on plots of Fi and its derivatives) that Fi{f3) is a com- 
pletely monotone function of P on (0,oo), i.e. {—!)'' d'Fi(P) / dp'' > for all /9 > 
and all integers k > O}^ Indeed, it seems that Gi(/3) = -Pi(/3) — 4//3 is completely 
monotone, which is stronger. 

Even more strongly, we conjecture (based on computations for Im/3 > 0) that 
G\{(3) = Fx{P) — 4//3 is a Stieltjes function for A = and A = 1, i.e. it can be written 
in the form 

/(/?) = ^ + / ^ (6-37) 

[0,oo) 

where C > and /) is a positive measure on [0, oo).^^ It is even possible that this 
holds also for < A < 1. □ 



Proof of Theorem 1.3. We modify the proof of Theorem 1.2 by using Proposi- 
tions 4.4 and 5.2 in place of Propositions 4.3 and 5.1, respectively. By Proposition 4.4 
[bound (4.10c)], for any x eV and S CV with a; G S*, we have 

\Cgis]{w)\ < TGis]{w)^{G[Slw)^\'\-'y' (6.38a) 

< TG[s]iw)'^iG,w)^^^\-^^/^ (6.38b) 

where w is defined by (4.12). We may apply Proposition 5.2 to deduce that 

\S\=n 



By the definition (4.12) of w we have 





\We\ 




\1 + We 


|l/2 



= A{G,w) (6.40a) 



A{G-x,w\g-x) < A{G,w') = A'{G,w) (6.40b) 



^^Scc e.g. [32] for the theory of completely monotone functions on (0, oo) in particular the 
Bernstein-Hausdorff-Widder theorem, which states that a function is completely monotone on (0, oo) 
if and only if it is the Laplace transform of a positive measure supported on [0, oo). 

^^More information on Stieltjes functions can be found in [32] [1, pp. 126-128] [2, 3, 4, 30] and the 
references cited therein. In order to test numerically the Stieltjes property for G'o(/3) and G'i(/3), we 
have used the complex-analysis characterization: a function / : (0, oo) ^ R is Stieltjes if and only if 
it is the restriction to (0, oo) of an analytic function on the cut plane C\ (— oo, 0] satisfying f{z) > 
for 2; > and Imf{z) < for Im^ > 0. See e.g. [1, p. 127] or [3]. 
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Hence 



S3X 
\S\=n 



. A(G,„r'Il±(^^ (6.41) 

where A = A'(G, w?) /A(G, «?). Therefore, the condition (3.4) for the weights (2.3) /(6.1) 
is verified as soon as 

inf (e" - 1)-^ M'^'^iG: w)^/^A{G, w)^-^ ^^^^(^Z ' - ^- (^-^2) 

n>2 ^' 



If we set L = \q\"^(G,w)~^^'^A(G,w)~^, this is precisely the inequahty contained 
in the right-hand side of (1.10a), or equivalently in the right-hand side of (6.5) 
with P — 1. So Zciq^w) whenever L > = Fx{l), i.e. whenever \q\ > 
Kl^{G,wy/^A{G,w). The equivalence with (1.10b) and the inequality (1.10c) 
then follow from Lemma 6.1. □ 

Remark. A weaker version of Theorem 1.3, in which and A(G, w) are replaced 
by Ki and A{G,w), respectively, can be proven by the same argument but using 
Proposition 5.1 in place of Proposition 5.2; moreover, we do not need to use (4.10c) 
but only (4.10d). □ 

Discussion. 1. We can now understand why the apparently minor improve- 
ment from Proposition 4.3 to Proposition 4.4 — changing a factor ^(G, ty)"/^ to 
^(G, to)^"~^^/^ — leads to the significant improvement (in most cases) of the final 
bound from Theorem 1.2 to Theorem 1.3, namely, replacing a growth ~ ^(G,iu) by 
~ \E'(G, wY^"^. Indeed, we can see using Lemma 6.1 that whenever we have a bound 
of the form 

^ \Cg[s]{w)\ < vl/(-i)/2+^' , (6.43) 



\S\=n 



we will obtain a bound on roots of Zciq, w) of the form 

|g| < A*^/^Fi(^-'') . (6.44) 

Proposition 4.3 and Theorem 1.2 correspond to 6 = 1/2, while Proposition 4.4 and 
Theorem 1.3 correspond to 6 = 0. 

2. Let us compare the bounds provided by Theorems 1.2 and 1.3: 

Theorem 1.2: IC* {G , w)) A' {G , w) (6.45a) 

Theorem 1.3: *(G', wf^^ A{G, w) (6.45b) 
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where A — A'(G, w)/A.{G, w) e (0, 1]. Their ratio is therefore 



Theorem 1.3 ^ Kl-^{G,wY'^ ^ F^jl) 

Theorem 1.2 A /C*(*(G', tu)) AFi(*(G,u?)-V2) " > 

Now, it is not difficult to see that A(G', w) < A'(G', w) ^(G, it))-*^/^, or in other words 
li?)^^/^ < A.^^ Since Fi[(3) is a decreasing function of (3 by Lemma 6.1(a), we 
have Fi(*(G','i/))-i/2) > Fi(A) and hence 

Theorem 1.3 Fxil) _ , . , . 

Theorem 1.2 ^ a7i(A) = '^^^ " ^'-''^ 

Both -Fa(I) and A-Fi(A) are increasing functions of A by Lemma 6.1(a,b), but their 
ratio g{X) does not have any obvious monotonicity. Numerically we find that g{X) 
decreases from the value Kq/4 pa 1.223222 at A = to a minimum value ~ 0.930714 
at A 3.70249, and then increases to 1 as A — > oo. We have not succeeded in proving 
that g{X) < g{0) for A e [0, 1], but if is true we can conclude that Theorem 1.3 is 
never more than a factor 1.223222 worse than Theorem 1.2. And in any case we 
have 

'^^^ - iJxZx) = X " ^-^26913 for A G [0, 1] . (6.48) 

We shall see in Examples 7.1 and 7.2 that Theorem 1.3 can indeed be up to a factor 
f» 1.223222 worse than Theorem 1.2. 

3. It is curious that the bound of Theorem 1.3 is not always better than that 
of Theorem 1.2, despite using better "ingredients" in its proof: namely, the bound 
(4.10c) from Proposition 4.4 always beats the bound (4.5b) from Proposition 4.3; and 
Proposition 5.2 is always better than Proposition 5.1. How is it that the final result 
can sometimes be worse? 

The explanation is that the ratio of the bounds (4.10c) and (4.5b), 

(410c) ^ T.(^)vI.(ff,^)(|vM)/^ ^ 

(4.5b) Th{w')^{H,w)\'^\/^ Th{w')^ ' ' ' ^ ' 

is the product of a "bad" factor Th{w)/Th{w') and a "good" factor ^!{H,w)~^/'^ . 
Now. the "l)ad" factor Tiiiw) /Th{w') is always bounded by ^(i7, lo)"^/^ — which is 



^^Proof. For each edge e we have 

mm.<\we\,-—^ PTTo > = min <^ We , -pT"^ r ^ x maxjl, 1 + } 



< min<|KI, ^r^^i> x ^{Cwf^ . 
Summing this over e3 x and taking the maximum over x €:V,we obtain the desired inequality. 
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why (4.10c) is always better than (4.5b) — so it follows that 



E Tgis]{w) 



S3x,\S\=n 



E Tgis]{w') 

S3x,\S\=n 



- < ^{G,wy^^ . (6.50) 



But there is no guarantee that the upper bounds on the numerator and denominator 
of (6.50), obtained by applying Propositions 5.2 and 5.1, respectively [cf. (6.39) and 
(6.3)], will also have a ratio < "^{G,wy^^. It is true that Proposition 5.2 is always 
better than Proposition 5.1 for a given set of weights w, but it could well happen 
that Proposition 5.2 applied to the weights w loses more than Proposition 5.1 applied 
to the weights w'. And this is indeed what happens in some cases (see Examples 7.1 
and 7.2). 

It is, nevertheless, somewhat disconcerting that Theorem 1.3 is not always better 
than Theorem 1.2. It would be interesting to try to find a single natural bound that 
simultaneously improves Theorems 1.2 and 1.3. 

4. Let us sketch how to obtain a bound that "interpolates" between Proposi- 
tions 4.3 and 4.4, and thence between Theorems 1.2 and 1.3, indexed by a parameter 
a G [0, 1]. The first step is to notice that the proofs of Propositions 4.3 and 4.4 bound 
the summand of (4.9b), individually for each tree T, in two different ways: one is 
(4.9c), and the other is the analogous bound in which w'^ = min{\we\, |we|/|l + 
is replaced by w" and the product is restricted to E \ E{x). So, taking for each T the 
weighted geometric mean of these two bounds, we obtain 

\Ch{w)\ < TH{w^''^)^{H,w)^\''\-''y^ (6.51) 

where 

fminjlwel, |t(7el/|l + "ifel^""^^} if 6 G E(x) 
[min{|i(;e|, |^ye|/|l + We|} otherwise 
Imitating the proof of Theorem 1.3, we then obtain the bound 

|g| < K{G,w)<i'{G,wY/'F,S'b{G,w)-^'-^y') (6.53) 

where 



K{G,w) = A(G,«;(»)) = max^minjl^.l, 'j^' | (6-54) 

e3x ^ I ^1 



(so that A(, = A' and A; = A) and A„ = A'{G, w)/A'^{G, w) G (0, 1]. We then have 
bound (6.53) Fa„(^(G, 'u;)-^^-")/^) 



Theorem 1.2 A„ Fi{^{G, wy^^) 



(6.55a) 



- *(G,M?)-«/2Fi(*(G,w?)-V2) ^ • ^ 
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where the inequahty uses \a > ^(G, it>)^"/^ together wirh the fact that F\{(3)/\ is 
a decreasing function of A by Lemma 6.1(c). Now, the upper hound (6.55b) is an 
increasing function of a e [0, 1] for all ^(G, w) >1: this follows from Lemma 6.1(c) 
with A = ^{G, II?) and ji — ^{G, w)"^/'^. But it is not clear what are the possible 
behaviors of the true bound (6.55a). For G — (see Example 7.1 below), wc have 
Aa = ^(G, lu)^"/^, so that (6.55b) is an equality, and the best bound is a = (i.e. 
Theorem 1.2). On the other hand, in cases like Examples 7.5 and 7.6, the best bound 
is clearly a — 1 (i.e. Theorem 1.3) because it has the smallest power of ^(G, w). We 
do not know whether there are examples in which some value a e (0, 1) might be 
optimal. □ 



Finally, let us prove Lemma 1.4 concerning the behavior of ^'(G, w) and A'(G, w) 
under parallel reduction: 

Proof of Lemma 1.4. Inequality (1.11) follows immediately from the fact that 
(1 + wi){l + W2) = 1 + ^^3. To prove (1-12), let us consider the following cases: 

Case 1 : ll+iwij < 1 and 11+^21 < 1. Then min ji^^j = I^Wil for 1 < i < 3, 

so we just have to prove that \w^\ < \w\\-\-\w2\- Since wz — w\-\-W2-\- W1W2, we have 

{Wsl = \Wi + 'W2 + WiW2\ = \Wi + 'W2{1 + Wi)\ < | Wi | + | u;2( 1 + ""^l) | 

= + 1^21 |1 + wil < Iwil + |'u;2| (6.56) 



since 11 + Wil < 1. 



Case 2: \1 + Wi\ > 1 and II + W2I > 1- Then min||wi|, j^^'^ 



II , I for 

1 < i < 3. Let w' = for 1 < i < 3, so that 1 + = (1 + Wi)'^ for 1 < i < 3 

— — i l-\-Wi — — ' % \ ^/ — — 

and hence (1 + w[){l + W2) = 1 + w'^. Since |1 + -u;^! < 1 and |1 + < 1, we may 
apply Case 1 to w[,w'2,w'^ to deduce that \w'^\ < \w'i \ + \w'2\, as required. 

Case 3: \1+Wi\ < 1, |l+ty2| > 1 and |l+t(;i| II+W2I < 1. Then min | j^^} = 

\wi\ for i e {1, 3}, and min ||w2|, |Y^} = By hypothesis we have |1 + t(;i| < 

|1 + 1(;2|~''^- Hence 

\W2 1 

\w3\ = l^i + ^2(1 + wi)| < |wi| + 1^2! |1 + wi| < |wi| + T- 7, (6.57) 

1 1 + W2 1 

as required. 

Case 4: \1+Wi\ < 1, |l+it;2| > 1 and \1+Wi\ \1+W2\ > 1. Then min | j^^} = 

\w,\, and min {1^,1, ^} = ^ for i e {2,3}. Let w', = for 1 < i < 3. 

Then jl + ty^l > 1 and |l + t(;2| < 1 with jl + ty^l jl + Wg] < 1, so we may apply Case 3 
(with indices 1 and 2 interchanged) to deduce that li^sl < jx^^ + |'"^2l — + l'"^2L 
as required. □ 
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Remcirk. We suspect that the transformation 



w' = --^ (6.58) 

1 + w ^ ^ 

employed in Cases 2 and 4, which satisfies (1 + w') = (1 + w)~^ and hence preserves 
the parallel-connection law (1 -\-Wi){l +'^^'2) = 1 + ^^3, may have other applications in 
the study of the multivariate Tutte polynomial. This transformation is involutive (i.e. 
w" = w), maps the complex antiferromagnetic regime |1 + w| < 1 onto the complex 
ferromagnetic regime > 1 and vice versa, and maps the real antiferromagnetic 

regime —l<w<0 onto the real ferromagnetic regime < w' < +00 and vice versa. 
In the physicists' notation w = e.'^ — 1 where J is the Potts-model coupling, the 
transformation (6.58) takes the simple form J' = —J, which makes its properties 
obvious. 

This fine of reasoning also suggests that the quantity 

A'(G,w) — max^^ min I Iwel, — '- — ri = max min{ | We | , (6.59) 

arising in Theorem 1.2, which we introduced simply because it arose naturally in our 
proof of Proposition 4.3, may also be "natural" in some more fundamental sense. □ 



7 Examples 

In this section we examine some examples that shed light on the extent to which 
Theorems 1.2 and 1.3 are sharp or non-sharp. For each graph G, we attempt to 
compute or estimate the quantity 

Qmax(G,io) = max{|g|: ^^(g, lo) = 0} (7.1) 

and compare it to the upper bounds given by Theorem 1.2 and Theorem 1.3. In what 
follows we abbreviate A'(G, to), 'kiG.w), '^{G,w), QmaAG,w) by A', A, Q^ax- 

Example 7.1 Let G = K2, where the single edge has weight w. Then Zk2{(1,'w) = 
q{q + w), so that Qmax = On the other hand, if |l + w| > 1 wc have A' = |t(7|/|l + 
w|, A = \w\/\l + w\^/'^, ^ = |1 + and A = A'/A = l/\l + w\'^^^. Theorem 1.2 gives 
the bound |g| < /C*(\l') A', which behaves like A\w\ as \w\ 00, while Theorem 1.3 
gives the bound \q\ < Kl'if^/^ A, which behaves like Kq\w\ a; 4.892888|w| as \w\ — >■ 
00. So Theorem 1.2 is off by a factor of 4 from the truth, while Theorem 1.3 is off 
by a factor of ^ 4.892888 from the truth. In particular. Theorem 1.3 is worse than 
Theorem 1.2 by a factor tending to K^/i ^ 1.223222. 
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For the special case of G = K2, the convergence conditions (6.4) and (6.42), which 
were used in the proofs of Theorems 1.2 and 1.3, respectively, become 

inf(e"-l)-V"|g|-^A'(G',M?)*(G',u?) < 1 (7.2) 

a>0 

inf(e"-l)'^e'"|gr^A(G,w;)*(G,w;)^/2 < 1 (7.3) 

a>0 

because the only polymer in the graph K2 has size n — 2. Since A'(G, w) ^(G, w) ~ 
A{G, w) ^'(G, wY^^ = \w\, we have 

(7.2) <^ (7.3) <^ \q\ > 4\w\ , (7.4) 

which differs from the truth Qmax = I'^^'l by a factor of 4. We can understand this 

behavior as follows: 

1) The lost factor of 4 comes from the fact that, for a polymer gas consisting 
of a single polymer S of cardinality l^*! = 2, the Gruber-Kunz-Fernandez-Procacci 
condition (Proposition 3.1) gives 5 7^ whenever \p{S)\ < 1/4:, whereas the truth is 
that S 7^ whenever \p{S)\ < 1. 

2) Though the convergence condition (6.4) involves a sum X]^2' terms for 
n > 2 make a negligible contribution in the limit |w| ^ 00 because with |g| of order 
Iwl we have 

|g|-("-^)A'(G',w;)"-^ *((?,«;)"/" = (|«;|/|g|)"-^ |1 + ^ (7.5) 

as \w\ — > 00 whenever n > 2. That is why Theorem 1.2 is off from the truth by the 
same factor 4 that we see in (7.4), despite the fact that its proof allows for arbitrarily 
large polymers that do not occur when G — K2. 

3) By contrast, in the convergence condition (6.42), the terms with n > 2 do not 
disappear in the limit |w| ^ 00 with |g| of order because 

[\q\-'A{G,w)^l>{G,wY/'r-' = i\w\/\q\r-' (7.6) 

is of order 1 for all n. This is why Theorem 1.3 is off from the truth by more than 
the factor 4 that we see in (7.4); we lose an additional factor K^/i pa 1.223222 by 
allowing for nonexistent large polymers. □ 

Example 7.2 In any simple graph G with at least one edge, we can choose weights 
w such that Theorem 1.2 beats Theorem 1.3 by a factor arbitrarily close to K^/AKi 
1.223222. It suffices to take Wg — w (with |1 + ^1 > 1) on all the edges of a nonempty 
matching, and We = wq on all other edges; then as wq ^ we have A' — >■ |t(;|/|l + wj, 
A ^ \w\/\l + w\^/^, ^ ^ |1 + w| and A = A'/A ^ l/\l+w\^/^. So the comparison 
of the bounds is the same as for G = K2, and Theorem 1.2 beats Theorem 1.3 by a 
factor tending to Kq/A ^ 1.223222 as \w\ 00. 
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For instance, let G be the n-cycle C„ with n > 3, taking We = w for exactly 
one edge and We = Wq for all other edges. Then Zo^q^w) = {q + w){q + Wo)"'~^ + 
wwQ^^{q — l). As |w| oo at fixed n and u^o, we have (5max(G, lu) = |w| +o(|w|). On 
the other hand, if jl + wl, jl + Wol > 1 we have A'(G', w) — \w\/\l + w\ + \wo\/\l + Wo\, 
A{G,w) = \w\/\l + w\'^/'^ + \wo\/\l + wo\^^'^ and *(G',ir) = \l + wo\ \l + w\. Therefore, 
as \w\ oo the bounds of Theorems 1.2 and 1.3 are 4(|wo| + |1 + '^^oDI'^^l + 0(1) 
and K^\l +wo\^/^\w\ + 0{\w\^/^), respectively, where ^ 4.892888. Both of these 
bounds have the correct magnitude as \w\ ^ oo ai fixed n and Wq. The bound given 
by Theorem 1.2 is better than that given by Theorem 1.3 when l^ol is small, and 
worse when l^ol is large. □ 

Example 7.3 Let G = ^ (a pair of vertices connected by k parallel edges) with 

We = w for all e. Then Za{q,w) = q[q+ (l + w)^ -1], so Qniax{G, w) = |(1 +-U7)^ - 1|. 
Now, if |1 + "U^l > 1 we have A'(G,w) = k\w\/\l + w\ and ^(G,w) = |1 + 
Therefore, as \w\ ^ oo at fixed k, the bound of Theorem 1.2 is a factor 4k from being 
sharp. 

On the other hand, we may first apply parallel reduction to yield a simple graph 
G — K2 with weight w — (l + w)'' — 1 on its single edge, and then apply Theorem 1.2 
or 1.3 to (G, w). The resulting bound is then (as \w\ 00) a factor 4 or 4.892888 
from being sharp (see Example 7.1). □ 

Example 7.4 Let G be the n-cycle C„ (which is simple for n > 3), with We = w 
for all e. Then Zdq.w) = {q + w)" + {q — l)w^. As \w\ ^ cxo at fixed n, we 
have Qma.yi{G,w) = IujI"/*^""-^) + 0(|w|). On the other hand, if |1 + w\ > 1 we 
have A'{G, w) = 2\w\/\l + w\, A{G, w) = 2\w\/\l + w\^/'^ and *(G, w) = \1 + wp. 
Therefore, as \w\ 00 the bounds of Theorems 1.2 and 1.3 are 81^1^ + 0(|u»|) and 
2Kq\w\^^^ + 0{\w\), respectively (here 2Kq ^ 9.785776). Both of these bounds have 
the wrong order of magnitude as j-u;! ^ 00 at fixed n > 4, but the bound given by 
Theorem 1.3 is a significant improvement over that given by Theorem 1.2. □ 

Example 7.5 Let G be the complete graph Kn- Take We = w > ior all e, with w 
fixed independent of n (unlike the usual [8] scaling w = X/n). Then Janson [18] has 
very recently proven that 

lim ^ log ^K„(e°", w) = maxfi log(l + w),a] for a > . (7.7) 

n—^oo n 

[This is because the sum (1.1) is dominated by two contributions: the terms with 
{V, A) connected, which together contribute e""(l + w)('^^[l + o(l)], and the term 
A — 0, which contributes e'*" .] It then follows from the Yang-Lee [35] theory of phase 
transitions (see e.g. [28, Theorem 3.1]) that Zx„(e""', w) must have complex roots q;„ 
that converge to = | log(l + w) as n ^ 00. Hence (5max(-f^n, w) > (1 + 'u;)"/^+''(") 
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(and this is presumably the actual order of magnitude). On the other hand, we have 
A{Kn,w) = (n - l)w/(l + and = (1 + so that the upper 

bound given by Theorem 1.3 is nearly sharp [it exceeds the truth by at most a factor 
0{n) even though both the truth and the bound are growing exponentially in n]. By 
contrast, the bound of Theorem 1.2 is much worse, because of its growth as ^{G, w) 
rather than ^(G, mj)^/^. □ 

Example 7.6 Let G be a large finite piece of the simple hypercubic lattice (for 
some fixed d > 2) with nearest-neighbor edges, and take We = w > Q for all e (here 
w > corresponds to the ferromagnetic case). For real g > sufficiently large, it is 
known [23, 22, 19, 21, 10] that the first-order phase-transition point Wt lies at 

wt{q) = q^'" + 0(1) . (7.8) 

It then follows from the Yang-Lee [35] theory of phase transitions that there will be 
complex zeros of the partition function arbitrarily close (as G grows) to the phase- 
transition point {q,Wt{q))\ so as w | oo (for fixed > 2) we will have asymptotically 
Qmax(G, w) > w'^[l + 0{l/w)] (and this is presumably the actual order of magnitude). 
Since A(G, w) = 2dw/{l +w)^^'^ and ^{G, w) = (1 +wy^, the upper bound given by 
Theorem 1.3 is off by at most a factor of order w^/^ (i.e. it gives a growth ~ ty'^+^Z^ 
instead of w'^). It is perhaps worth observing that the bound would be off by only 
a bounded factor if Theorem 1.3 could be improved to use A'(G, lo) rather than 
A(G, w). By contrast, the bound of Theorem 1.2 is again much worse, because of its 
growth as ^(G, w) rather than ^(G, ti;)^/^. □ 
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